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We present experimental results on the optical absorption spectra of epitaxial graphene from the 
visible to the terahertz (THz) frequency range. In the THz range, the absorption is dominated by 
intraband processes with a frequency dependence similar to the Drude model. In the near IR range, 
the absorption is due to interband processes and the measured optical conductivity is close to the 
theoretical value of e^ /Ah. We extract values for the carrier densities, the number of carbon atom 
layers, and the intraband scattering times from the measurements. 



Graphene is a single atomic layer of carbon atoms 
forming a honeycomb crystal lattice P, ^ . The unusual 
electronic and optical properties of graphene have gen- 
erated interest in both experimental and theoretical are- 
nas 0, 13, 0, S S|- The high mobility of electrons in 
graphene has prompted a large number of investigations 
into graphene based high speed electronic devices, such 
as field-effect transistors and pn-junction diodes, pho- 
tonic devices, such as terahertz oscillators, and also low 
noise electronic and optical sensors [1, 0, S, 0, [13, U | • 
For many of these applications, knowledge of the optical 
properties of graphene is critical. 

Graphene layers can be obtained via micromechanical 
cleaving (exfoliation) of bulk graphite followed by careful 
selection of monolayers by using optical, atomic force, or 
scanning electron microscopes Although this tech- 
nique results in relatively high quality films, it might not 
be suitable for large scale production. Recently, epitaxial 
growth of graphene by thermal decomposition of SiC sur- 
face at high temperatures has been demonstrated [l, • 
This technique can provide anywhere from a few mono- 
layers of graphene to several (> 50) layers on the surface 
of a SiC wafer. Graphene layers grown by this technique 
have demonstrated low temperature carrier mobilities in 
the few tens of thousand cm^/V-s range [6J. In addition, 
the electronic as well phononic properties of epitaxially 
grown graphene multilayers have been found to be differ- 
ent from those of bulk graphite and similar to those of 
a graphene monolayer [ 



13, 14, 15, n 



This observed 
difference in the properties of epitaxial graphene and 
bulk graphite has been attributed to a different stack- 
ing scheme for carbon atom layers in epitaxial graphene 
compared to that in bulk graphite in which the layers are 
stacked according to the Bernal scheme 3, \mM- Al- 
though the band energy dispersion in carbon atom multi- 
layers can exhibit massless Dirac Fermion like behavior at 
small energies for various different stacking schemes and 
interlayer couplings, the band energy dispersions at large 
energies can be significantly different [l3, [l3| • The exact 



structure of epitaxial graphene and the nature of inter- 
layer couplings remain active areas of investigation. Mea- 
surement of the optical absorption spectra over a wide 
frequency range can provide useful information about the 
structure of epitaxial graphene. 

In the visible to the mid-IR wavelength range (A < 
10 /im), the optical absorption spectra of exfoliated 
graphene monolayers has been reported recently [13, [H, 
,19]. In this paper, we report results from measurements 
of the optical absorption spectra of epitaxial graphene 
from the visible to the terahertz frequency range for the 
first time and compare the results with the theoretical 
predictions for graphene. In graphene, the valence and 
conduction bands resulting from the mixing of the p^- 
orbitals are degenerate at the K (K') points of the Bril- 
louin zone. Near these points, the tight-binding Hamil- 
tonian in the nearest-neighbor approximation , using the 
basis consisting of orbitals centered on the A and B 
atoms, can be linearized and written as [l[. 



H = 



A hv [kx + i ky) 

h V {kx — i ky) —A 



(1) 



where v « 10^ m/s is the Fermi velocity. This 
Hamiltonian results in the energy dispersion relation for 
the conduction and valence bands given by, Ec{k) = 



±-\/A2 ~\- {Tiv ky. The bandgap is equal to 2 A and 
could acquire a non-zero value as a result of any inter- 
action that breaks the symmetry between the A and B 
atoms in the unit cell of graphene. Optical absorption in 
graphene is described by the optical conductivity (t{lu). 
It can be written as the sum of the interband conductiv- 
ity cTinterii^) and the intraband conductivity Ointrai'^), 
both of which can be found using the Hamiltonian above 

13, [13. 



and are given below 
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Here, /(e — Ej) is the Fermi distribution function with 
Fermi enery F describes the broadening of the in- 
terband transitions, and r is the momentum relaxation 
time due to carrier intraband scattering. The frequency 
dependencies of the real parts of Ointeriy) and <7intra{i^) 
are depicted in Fig.l, assuming A = 0, F = 10 meV 
and T=300K. Fig. 1(a) shows the conductivities for Ef = 
— 100 meV and two different values of the scattering time 
r: 25 fs and 5 fs. Fig. 1(b) shows the conductivities for 
T equal to 25 fs and two different values of the fermi 
energy Ef. meV and -100 meV. At large frequencies, 
the real part of tTmter(w) has a constant value equal to 
e^/Ah. At small frequencies, the real part of (Tmter(w) 
approaches zero because interband optical transitions are 
blocked due to the presence of electrons and holes near 
the band edges. The plasmon dispersion and the free- 
carrier absorption in graphene are described by Ointrai'^)- 
Its frequency dependence is similar to that of a Drude 
model, as it is evident from the pre- factor in Eq. ^ . The 
real part of (Jintra{i^) approaches zero for large frequen- 
cies at which carriers are unable to respond. At small 
frequencies, the real part of (Ti„tra(w) approaches the DC 
conductivity of graphene. Fig.l (b) shows that the spec- 
tral shape of crintrai'^) at small frequencies is strongly in- 
fluenced by the intraband carrier scattering time r. The 
total conductivity <y{i^) has a minimum in the frequency 
range where both crintrai'^) and amteri^) arc small. 

The epitaxial graphene samples used in this work 
were grown on the carbon face of semi-insulating 6H- 
SiC wafers using the techniques that have been re- 
ported in detail previously [l^ . The samples were grown 
at temperatures of 1400 °C - 1600 °C and pressures of 
2 — 7 X 10~^ torr. The number of carbon atom layers in 
each sample were estimated through X-ray photoelectron 
spectroscopy (XPS) using the Thickogram method 
Raman spectroscopy (using excitation wavelength of 488 
nm) of the samples showed a single-resonant G peak close 
to 1580 cm~^, a double-resonant D' peak close to 2700 
cm^^, and also a relatively low intensity double-resonant 
D peak near 1350 cm~^ 0, HH. The D peak is not al- 
lowed in perfect graphene layers since it requires an elas- 
tic scattering process, which is made possible by disorder, 
to satisfy momentum conservation [13| . The presence of 
the D peak, therefore, indicates the presence of disorder 
in the samples. The ratio of the intensities of the G and 
D peaks [Ig/Id) has been shown to be proportional to 
the crystal coherence length [2^. The ratio Ig/Id for 
samples A, B and C were ~13, ~17 and~2 respectively, 
indicating that sample C has a much larger level of disor- 
der compared to the other two samples. Ultrafast carrier 
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FIG. 1: Real parts of the interband (solid) and intraband 
(dashed) optical conductivities of graphene are plotted, (a) 
Ef = —100 meV and r equals 25 fs and 5 fs. (b) r = 25 fs 
and Ef equals meV and -100 meV. Values of F and A are 
assumed to be 10 meV and 0, respectively, and T=300K. 



dynamics in samples A and B have been studied in a pre- 
vious work ^T^. Sample C reported in this paper is not 
the same as the one reported in our earlier study 22 1. 

Three different instruments were used to measure the 
optical transmission through the graphene samples. In 
the visible to the near-IR wavelength range (0.4-0.9 /xm) 
a grating spectrometer was used. In the near-IR to the 
mid-IR range (1.4-25 /zm) a mid-IR Fourier Transform 
IR (FTIR) spectrometer was used. And in the mid-IR 
to the far-IR (terahertz) wavelength range (15-200 /im) 
a far-IR FTIR spectrometer was used. The measured 
transmission spectrum for each sample was normalized 
to the transmission spectrum of a reference SiC wafer. A 
small (few mm size) aperture was used to ensure trans- 
mission through equal areas of the sample and the ref- 
erence and also to avoid corruption of data when chang- 
ing instruments due to possible non-uniformities in the 
sample. The SiC substrate transmits very little in the 
6-14 jxm wavelength range due to multi-phonon absorp- 
tion 24'! . As a result, the measured transmission spectra 
had poor signal-to-noise ratios in this wavelength range. 
The fringes in the transmission spectra arising from mul- 
tiple reflections within the SiC substrate were numeri- 
cally filtered out after normalizing with respect to the 
transmission spectrum of the reference SiC wafer. 

Fig. 2 shows the normalized transmission spectra 
through the epitaxial graphene samples A, B and C. The 



3 




Wavelength (|xm) 




N = 12 
T = 20 fs 
Ef = - 95 meV 



1 10 100 

Wavelength (|a.m) 




Sample C; 

N = 20 
x = 4 fs 
Ef = - 75 meV 



1 10 100 

Wavelength dim) 



FIG. 2; Measured transmission spectra (solid gray lines) of 
samples A, B, and C from 0.4 fim to 200 /jm along with 
the theoretical transmission spectra (dashed black line) using 
Eqs.([2}-(|4]). The values of the fitting parameters are included 
in the insets. The value of A is assumed to be zero. 



optical absorption in the graphene layers depends sen- 
sitively on the substrate index of refraction. Matching 
the optical boundary conditions at the air/graphene/SiC 
interfaces, the optical transmission T(uj) through N 
graphene layers on a SiC wafer (normalized to the trans- 
mission through a plain SiC wafer) can be written in 
terms of the complex optical conductivity cr(a;) as, 



T{u;) = |1 + Na{Lu)^fio/eo/{l + ns^c)\' 



(4) 



where nsic ^ 2.55 is the refractive index of SiC. For 
N\(T{uj)\y^ fj.o/^o/i'i^ + nsic) « ^,T{u;) is related to only 
the real part of the optical conductivity and can therefore 
be used to measure the absorption spectra of graphene. 
We have used Eqs.([2|)-(|4]) to model the measured trans- 
mission data using N, t, Ef as the fitting parameters. 
The parameters F and A were assumed to be '-^ 10 meV 
and ~ meV respectively. Changing the values of F 
and A by small amounts (much less than ksT) had lit- 
tle effect on the final results. The dashed black lines in 
Fig. 2 are the theoretical fits to the experimental data 
(solid gray lines). The values of the fitting parameters 
are shown in the insets. A 5-10% variation in any one 
of the parameters, N, r,and Ef, produces a noticeable 
degradation in the quality of the theoretical fit to the 
data. 

The number of graphene layers N obtained this way 
agrees well with the value obtained through XPS. For 
example, the XPS method gave values of N equal to 6 
and 11 for samples A and B respectively. The extracted 
values of the Fermi energy correspond to average carrier 
densities of 5 x 10^^ cm~^ and ^ 8x10^^ cm~^ per layer 
for samples A and B, respectively. Due to the electron- 
hole symmetry of the graphene bandstructure near the 
band edge, both negative and positive signs of the Fermi 
energy will fit the experimental data equally well. Also 
note that the experiments only give information on the 
total conductivity Na{uj) of all the layers, and therefore 
the extracted value of the Fermi energy should be taken 
as an average value for all the layers. Recent work on epi- 
taxial graphene suggests that a concentration of carriers 
larger than the intrinsic value is expected only in the first 
few carbon atom layers near the SiC interface. Assuming 
that only the first two layers have non-zero Fermi levels 
and the remaining layers are intrinsic, values of the Fermi 
level equal to -150 meV and -290 meV for the first two 
layers of samples A and B, respectively, provide good fits 
to the measured data. As mentioned earlier. Sample C 
is significantly more disordered than Samples A and B. 
The transmission spectra of Sample C shows a distinctly 
different shape in the THz region (Fig. 2c), which can be 
fitted well with a very short carrier scattering time of ~4 
fs. 

It has been recently pointed out that the first few 
graphene layers in epitaxially grown graphene could ac- 
quire a bandgap as a result of interaction with the atoms 
in the SiC substrate [1^. Although a value of A equal 
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to zero was found to fit our measured data well, a value 
of A much smaller than fc^T would be difficult to detect 
in our measurements. A non-zero value of the bandgap 
2A in the mid-IR to far-IR range, where the intraband 
contribution to the conductivity dominates, would have 
little effect on the transmission spectra. If the value of 
the bandgap is in the near-IR to the mid-IR range, its 
effects at room temperature would be hard to distinguish 
from the reduction in the interband conductivity at small 
frequencies due to band filling effects (see Fig. 1(b)). Also, 
in multilayer graphene structures the optical response is 
dominated by the large number of layers that are not 
close to the substrate and do not have a bandgap. 

The short wavelength end of the measured transmis- 
sion spectra in Fig. 2 deviates from the theoretical pre- 
dictions for wavelengths shorter than ~ 2.5 fxm. The 
deviation is minimum for sample B and corresponds to 
^50% more absorption at 0.4 fim compared to the theory. 
The reasons for this deviation are not clear. Two factors 
could be responsible for this behavior: (i) the band en- 
ergy dispersion and interband optical matrix elements 
at large energies are different from those obtained from 
the Hamiltonian given in Eq.([2), and (ii) increased light 
scattering may be expected from the sample at shorter 
wavelengths as the wavelength approaches the crystal co- 
herence length (^^50-100 nm). Using a full-band tight- 
binding model for a graphene monolayer that includes 
second and third neighbor interactions [26| , and using the 
method for calculating the interband optical matrix ele- 
ments described by Johnson et. al. [13], we obtain values 
of the optical interband conductivity that deviate from 
e'^/Ah only for wavelengths shorter than 1.5 /im and are 
only 10% — 20% larger than e^/ih at wavelengths close to 
0.4 /zm. Therefore, effects related to trigonal warping and 
deviation of bands from linearity cannot be completely 
responsible for the observed discrepancy. The band en- 
ergy dispersion in epitaxial graphene at large energies 
could also be affected by the nature of the interlayer cou- 
plings. Note that the deviation of the measured trans- 
mission spectra from the theory is not the same for the 
three samples indicating that disorder might also have 
a role to play. More work is needed to investigate the 
nature of this discrepancy. 

In conclusion, we have measured the optical absorp- 
tion spectra of epitaxial graphene from the terahertz to 
the visible frequencies. The experimental results have 
been shown to be in agreement with the theory except 
at short wavelengths. Our results confirm the Drude-like 
frequency dependence of the intraband conductivity of 
graphene in the THz frequency range. The results pre- 
sented here indicate that absorption spectroscopy can be 
used as a noninvasive technique to characterize graphene 
films and find the values of parameters, such as the Fermi 
energy and the carrier density, carrier intraband scatter- 
ing time, and the number of graphene layers. The authors 
acknowledge support from the National Science Founda- 
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